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*^ ' Abstract 

(N; 

I In this paper, we show the integrabihty of spin- 1/2 XXZ Heisenberg chain with two 

arbitrary spin boundary Impurities. By using the fusion method, we generahze it to the 
spin-1 XXZ chain. Then the eigenvalues of Hamiltonians of these models are obtained by 
the means of Bethe ansatz method. 
^ ■ Pacs: 75.20H, 12.40E, 75.10D 

^-H ■ 1. Introduction 

■^ , Recently, more and more papers have focused on the Kondo problem. It is well known that 

Cd I the spin dynamics of the Kondo problem is equivalent to the dynamics of the spin chain with 

^ ■ magnetic impurities O. Although magnetic impurities play important role in the model, they 

F^ . usually destroy the integrability of the system. So how to maintain the integrability of the 

P5 I quantum impurity system is an important problem. Fruitful achievement has been obtained 

O ■ based on the methods of the Bosonization and renormalization technique, conformal field 

• • . theory and exact diagonalization in this field [2-5]. 

.!^ I The quantum inverse scattering method and the Bethe ansatz technique have been pow- 

p\ ' erful tools to study the integrable impurity problems in 1-dimensional physical system, such 

j^ ■ as Wang et al 's papers about the vertex model[6-8] and Frahm et al and Links et ai's series 

papers about the t — j model with impurities [9- 15]. The Heisenberg chain is an important 
model in the integrable system so many papers pay their attention to it. Andrei and Jo- 
hannesson first considered the integrable Heisenberg chain with impurities under periodic 
boundary condition |16|. Then Lee and Schlottmann generalized their results to arbitrary 



spin impurities [17, [I^. But they have to present some unphysical terms in the Hamiltonian 
to maintain its integrability, though those terms may be irrelevant [^ . To the open bound- 
ary condition problem, Gaudin considered the nonlinear Schrodinger model and the spin-1/2 
Heisenberg chain with simple open boundaries [ 19 1, then Schulz and Alcaraz|2C, ^]et al. gen- 
eralized it to Hubbard and other models. Wang have discussed the properties of the impurities 
with arbitrary spin coupled to the spin-1/2 XXX chain [B|. The spin-1/2 XXZ chain coupled 
with spin-1/2 impurities has also been discussed in Ref.[0]. In Ref.lg], the integrability of the 
spin-1 XXX chain with arbitrary spin impurities has been investigated. 

In this paper, we study the integrability of the open Heisenberg chain coupled with ar- 
bitrary spin impurities. We discuss the spin-1/2 XXZ chain in the first part of the present 
paper. The spin-1 case is presented in the second part. A brief discussion about our results 



is given in the last section. 

2. The spin-1/2 XXZ chain 



The i?-niatrix of spin-1/2 XXZ Heisenberg chain can be written as 

n Clin II din n fi 

R{u) 



( sin(n + 7?) 

sinti sinr/ 

sinr/ sinu 

sin(u + 7?) ) 



(1) 
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This i?-niatrix is regular and satisfies the unitarity condition 

R{u)R{—u) = sin(n + rj) sin(— n + f]) = p{u). 

If we suppose the first and the second space of the i?-matrix are auxiliary and quantum space 
respectively, and this i?-matrix, as an operator matrix , can also be written as L-operator 
form 



Ln{u) 






where 
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wa — w^ = sinu, 

W4, + w-i = s\n{u + rj), 
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The L-operator and i?-matrix satisfy the following Yang-Baxter relations (YBR) 



(2) 



Ri2{u - v) Li (u) Li (v) =Li (v) Li {u)Ri2{u - v), 

12 2 1 

Ruiu -v)T (u) T (v) =T (v) T {u)Ri2{u - v), 
where R12 acts on the auxiliary space f 1 t;2 , and the T defined by 

N 



Tiu) = l[L,< 



u 



1 



(3) 



acts on quantum spaces as wi ® t;2 ® • • • (X> f at. Here we have used the notation A= A0 1, 

A=l(^A. 

In order to construct the open boundary condition consistent with the integrability, we 
consider the refiection equation 

Ruiu -v) K {u)R2i{u + v) K {v) =K {v)Ri2{u + v) K {u)R2i{u - v). (4) 



where K is the reflecting matrix, which determines the boundary terms in the Hamiltonian. 
One can prove that the double-row monodromy matrix defined by U{u) = T{u)K{u)T~^{—u) 
also satisfies the reflection equation 

Ri2{u -v)U {u)R2i{u + v)U (v) =U {v)Ri2{u + v)U {u)R2i{u - v). (5) 



The dual iT-matrix K^[u) can be defined by the automorphism [ 26 1 

<j):K{u)^K+{u)=K\-u-r]). (6) 

It satisfies the dual reflection equation 

1 2 

Ri2{-u + v) K+ {u)Ri2{-u -v-2r]) K+ {v) 

2 1 

= K+ {v)Ri2{-u -v-2r]) K+ {u)Ri2{-u + v). (7) 

Then the transfer matrix can be defined by 

t{u) = trK+{u)U{u). (8) 

One can check that it satisfies the commutation relation 

[t(n),tW]=0. (9) 

Now we couple the spin-1/2 XXZ chain with two arbitrary spin impurities located at the 
ends of the system. Then the L-operator of the boundary cites can be written as 

, . / sin('u + S + df ??) d:[ s'mr] \ 

^^(")=( d+sin, sin(. ; f - df ,) j ' (^ = "' ^) (10) 

where d ,d^ are components of an arbitrary spin m of SUq{2). One can easily check that 
the L-operator satisfies the first relation of (3). It also have the unitarity relation 

Li{u)Li{-u) = d^ sin^ r/ -|- sin(u -|- -) sin(-ti + :t) = Pd{u) 

with d^ = sin(?Tir/) sin(r/ -|- mrj)/ sin^ rj. 
Define 

T{u) = Lb{u + Cb)LN{u)- ■■ L2{u)Li{u)La{u + Ca), 

f(u) = T-^{-u) X const. 

= La{u-Ca)Li{u)L2{u)---LM{u)Lb{u-Cb). (11) 



where Cj are free parameters. According to Cherednik[p5|| and Sklyanin's work pq| , the re- 
flection matrix and its dual are defined by 

x^/ X ,. /, sin(f — n)\ ^^,, , ,. / sin(f + -I- n -I- r?) \ ,^. 

K{u) = diag 1, . ;^ , K+{u) = diag 1, -^ ^ . 12 

V sm(^ + n)/ y sm(4+ — n — r/)y 

Recalling the definition of t{u) (8), one can check that the above formulas satisfy the com- 
mutation relation (9). By expanding t{u) in terms of u, we can obtain infinite number of 



conserved quantities which ensures the integrabihty of the model. The Hamiltonian of this 
model can be written as 

rr ^ 1 ^,dt{u) _ 

2p2^(0)p,(c„)pd(cb)tr,i^+(0) du '"-° 

N-l TT ( \ 

- Z^ IT— l«=o 

1 1 

, d{La [u^Cg) Liu-Cg)) . 1. . ^ dKx(u) \ 

^ 7, — i—T-j M=o+ L (n + Co)- — , . , L\u- Cg)\u=Q 

2pd{ca)du 2pd{cg)du 

'^ T T 

trr K+ (u) Lb {u + Cb)Hr,N{u) Lb (n - Cb) . 

H 7 T Im=o 

Pd{cb)p^{u)trr K+ (u) 

'^ T T 

trr K+ {u)d{Lb {u + Cb) Lb {u - Cb)) , , , /, o^ 

H 7 |n=o + const. (13) 

2pii{cb)trr K+ {u)du 

where Hjj^i{u) = — ^'^^^ Rjj^i{u). Denoting by Tj the ith term of the right hand side of 
(13), we have 

, N-l 

^1 = — : E (^1 • ^1+1 + ^1 • ^i+i + cos W^ ■ <^i+i 

T2 +n = ^{1 + af) ■ Ag + ^{1 - af) • 5„ + a+ • Cgdg + af • d+C^, 

r4 + Ts = ^(1 + (7^) • A + ^(1 - ^1/) -Bb + a^- Cbdl + (Tjv • d+a, 

1 /^„:„^_„ , o,,2„X 2 COS ^2 • 2, 



Ai = -— sin(?7 + 2d''ri) —id svcrirj) - sinid'^r]) s'mir] + d^r/)) 

2pd{ci) \ sm^ 

Bi = -— ( sin(77 - 2d''r]) ^^^ sin(ci + - - d'^'n) sin(-Cj + - - d'^r]) ] , 

Q = ^^^ sin(Q + ^ - e + d^v), (i = a, b) (14) 

where ^ should be changed to ^^ when i = b, and this Hamiltonian is hermitician when we 
choose pure imaginary q. 

To construct the algebraic Bethe ansatz, we rewrite the double-row monodromy matrix 
U(u) in the form 

Using the reflection equation (5), we can obtain the following commutation relation 

B{u)B{v) = B{v)B{u); 

A{u)B{v) = ^— —-^- -B{v)A{u) 

sm(7? + v + u) s\n{v - u) 

sin77sin(u + v) „, x a, x 

sin(77 + t; + n) sm(ti — v) 



sin(77 + t; + n) 



V{u)B{v) = ^- ^-^- r — B{v)V{u) 

sm{u + v + 7]) sm(n - v) 



sm(u — v) sin(2w + 77) 

sm7?sin(2t;)sin(2n + 2?7) 

+ ^ r^^ ^^WAiv), (16) 

sm(2i; + 7?) sin(u + t; + 77) 

where X'(u) = sin(2n + rj)T>{u) — sm.r]A{u). Using the relation (8) and (15), the transfer 
matrix t(u) can now be written as 

/ N +^/ X + 4/ X sin(f+ + u + 77) ~ , , sin(£+ — li) sin(2n + 277) ,, , , „, 

sm(2n + 77) sin(zti + r/) 

Define the pseudo- vacuum state |0) 

a+|0)=d+|0) = 0, (i = l, 2, •••, N) (18) 

Acting the elements of Ur{u) on |0), we have 

Ciu)\0) = 
A{u)\0) = wi\0) 

= sin(,^ + n) sin ("U + ^) _[| Jl sin(n + rcj + 7//2 + 77177) |0) 

r=±l j=a,6 

V{u)\0) = W2\0) 

= sin(^ — « — r/) sin(2u) sin^ ^ IT IT sin(ti + rcj + r//2 — mr/)|0) (19) 

r=±l i=a,b 

The eigenstates of t{u) can be constructed from the pseudo- vacuum state 

M 

\n) = llB{um- (20) 

Thus we obtain the eigenvalue of t{u) acting on |i7) 

,/ Mo\ +TT sin(77 + 7;i -u)sin(u + 7;i) 

1^-^ sm(r/ + -yj + m) sm(— n + Vj) 

+ t.,t.+ -fT ^^"^\-"- + ^^""^" + "- + '^V ) (21) 

rl^-": sin(« + f « + 77) sin(u — Uj) 

with the Bethe ansatz 

smjvj - g+ - ^) sin(7;j + ^ - ^) sin2^(7;j + 2) 



sin(t;, + e+ + 2) sin(7;,- - ^ + f) sin2^(7;,- - f, 



T-r T-r sin(-7;j + rck + ttit/) _ -j-j- sin(t;j - -Uj + r/) sin(7Jj + fj + 77) 
,i|i fcif^fe sin(t;j + rcfc - mr]) f-J-. sm{vj - Vi - r]) sm{vj + Vi - r])' 



3 The spin-1 Heisenberg chain 



3.1 The Zamolodchikov-Fateev 19-vertex model 



The Zamolodchikov-Fateev 19-vertex -R-matrix[p^] associated with the spin-1 representa- 



tion of ?7q(g/2) p3[| can be obtained by using the fusion method[22, 24]. It reads 
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(23) 



with 

ai = sin(u + 2r]) sin(n + 7]), a2 = sin u sin(n + rj) , 

as = sin(2r/) sin(ti + r/), a^ = sinusin(ii — rj), 

as = sinMsin(2r/), ag = sinr/sin(2?7), 

aj = a2 + ae- 
This i?-inatrix is regular and satisfies the unitarity relation 

sin(u -|- ry) sin(n — r/) sin(u -|- 2r/) sin(u — 2ry) = p{u). 

It satisfies the Yang Baxer equation 

Ruiu - v)Ri3iu)R23{v) = R23iv)Rniu)Ri2iu - v). (24) 

Inami et al^^ have obtained the general solution K{u) of (4). In this paper, we only adopt 
its diagonal form as in the Ref.pc 



K{u) 



diag{ki{u),k2{u),k3{u)) 

/ sin(§-|-^-n) sin(§ +^ - u)sin(2 - ^ + m) 
diag 1, 



sin(2 + ^ + n) ' sin(2 + ^ + u) sin(2 - ^ 
and the corresponding dual reflection matrix takes the form 
K+{u) = diag{k+{u),k^{u),k+{u)) 



(25) 



diag 1, 



sin(f + e+ + u) sin(f + e+ + n) sin(-2 - ^+ - u) 



sin(-2+eH 



u)'sin(-2+^H 



u)sm{^-^++u)^ 



(26) 



3.2 Fusion of the boundary L-operator 



In this section we discuss the fusion procedure of the boundary L-operator. The permu- 
tation operator and the projection operators are defined by 



Pi 



12 



P. 



12 



P 



12 



i?i2(0)/sin(r/) 

-Pi2(-r/)/(2sinr/) 

1-Pi 



12 



(27) 
(28) 
(29) 



respectively, where the R{u) is from the 6- vertex model (1). They satisfy the following 
properties 

A2 = 1' (-P12) — -^^12' 

P+Pf2 = Pf2P+ = 0. (30) 

Now we use fusion procedure to obtain the high-dimensional L-operator. Taking v = u + rj 
in the equation (3), the YBR gives 

Ri2{-ri)Lid{u)L2d{u + r]) = L2d{u + 7])Lid{u)Ri2{-ri), (31) 

i 

where d represent the boundary terms, and we write Lj as Lij for convenience. Multiplying 
above equation by P{2 from the left and right respectively, we get 

Pt2L2d{u + r])L^d{u)P{2 = 0, (32) 

P{^LUn)L2d{u + r,)P+ = 0. (33) 

Define 

L<i2>d{u) = Pt2Lid{u)L2d{u + r])P+, (34) 

L'<i2>diu) = P^2L2d{u)Lid{u - v)Pu^ (35) 

which satisfy the YBR respectively. 

^<12><34>(^t - v)L^i2>d{u)L^s4ydiv) = i><34>d(t')^<12>d(^i)^<12><34> (^t " v), (36) 

^<12><34>(^ - v)L'^i2>d(.u)L'^s4yd{v) = ^<34>d(^')^<12>dW^<12><34>(^ - v), (37) 

where i?<i2><34>(^ — v) is the fused i?-matrix|24], acts on l^<i2> '53 ^<34>- Here we give the 



proof for R(u): 

LHS = R<i2><3'i>iu-v)P^2Lidiu)L2d{u + r])P^2Pu^sd{v)L4d{v + 7])P;^4^ 
= P^^P^^Riiiu -V- r])Ri3{u - v)P^iPtiR2A{u - v)R23iu -v + r]) 

xP+P+Lid{u)L2d{u + i])Lm{v)Lm{v + ri)Pi2Pu 
= Ri4,{u -V - r])Ri3{u - v)R2a{u - v)R23{u -v + ri)Lid{u)L2d{u + r]) 

xL3d{v)Lid{v + v)Pi2Pu 
= Ri4,{u -V - ri)R2A{u - v)Ri3{u - v)Lid{u)L3d{v)L2d{u + rf) 

XR23{U -v + r])L4d{v + r?)-Pi2-f34 

= L3d{v)L4div + r])Lid{u)L2d{u + r])Ri4{u -v- r])Ri3{u - v) 

XR24{U - v)R23{u -V + v)PuPl2 

= PtiL3d{v)Lid{v + r])P+P+Lu{u)L2diu + 7j)P+P+P+Ru{u - v - r]) 
xRi3{u - v)P+4P+4R2i{u - v)R23{u -v + v)P3aPi2 

= L^3iyd{v)L<l2>d{u)R<12><M>{u - v) 

= RHS. (38) 

The proof for the other formula is similar. Substituting relation (10) into (33) and taking 
the transformation 



/ 1 



L<,i2>d{u) :^ (1 V2cosr/ l)L<i2>d(n) 



\/2 cos r] 

V 1 



^ ■ (39) 



we have 



with 



an 


ai2 


ai3 


021 


022 


023 


031 


032 


033 



L<i2>d{u) = a2i a22 023 , (40) 



an = sin{u + d^ rj) sm{u + 7] + d'^ 7]) , 

012 = d^ \J2 cos T] sin(n + d^T]) sin ry, 

013 = (d")^sin^r/, 

a2i = v^2 cos T] sin(u + d^?]) sin r\d^ , 

022 = sin u sin(n + r/) + d sin 77 — 2 sin (d^T]) cos 77, 

a23 = v^2 cos ry sin(u — d^r\) sin ryd" , 

031 = ((i+)^sin^?7, 

032 = d'^ \J2 cos 77 sin(n — d^rj) sin ry, 

033 = sin(n — d^rj) sm.{u + 7] — d^rfj. 

This L-operator also satisfies the unitarity relation with 

Pd{u) = t{(— rf — cos(2n — 7?) + COST? + d cos(2?7)} 
x{—d — cos{2u + r]) + cos r] + d cos{2ri)}. 



3.3 The Hamiltonian of this model 

Define 

T{u) = Lb{u + Cb)LNiu)---L2{u)Li{u)La{u + Ca), 

f{u) = T~^{-u) X const. (41) 

= La{u- Ca)Li{u)L2{u)---LNiu)Lb{u- Cb), (42) 

where Ca and Cb are constant. The spin-1 L-operator is obtained from the i?-matrix (23) by 
assigning the second space to be the quantum space, and Li{u) {i = a, b) is given by (40). 
The Hamiltonian of this model is as same as the spin-1/2 case (13). Here we give Tj as 

Ti = -J— J: I ^^^^ - ^'^ ' T^' + (1 - cos r?) (44+i4.-,i + s-s^, 
sm(277) -^ I COST] cos^ r/ \ j j^ j j-r^ j j^^ 

- (1 - cos(27?)) (s|.|+i - (sf)^(.|+i)' + (.|)' + {sUf)} 

I ^ ^ 1_ 1__ sin(s|?7) sin(sL 1??) 

Sj • Sj+i = -s- sj^-^ + -s- s-^^ + cosr/ '—^ — '- , 

COST? 2 ■' ■'^ 2 •' ■'^ sm ?7 

oio\ /ooo 



with 



and 



^J2cosr| I 1 I , s = y/2 cos r/ 




1 
\ 1 




To + T. 



-{s^ + si)Aa + s\ni]s^s+Dad + {s+YEa{d 



7-^2 



+ sin?7s Szd^Da + 



1 



2 cos Tj 



(s+s +s s^) - s^^Ba-smrjs'^SzFad 



+ {s-Y{d+YEa - smr]s,s'd+Fa + -{si - S,)Ca, 



T4 + T5 



1 



{s^ + sl)Ab + sin r]s^s+Dhd + {s+fEbid' 



with 



Bi 



C, 



Di 



Fi 



+ sin 77s Szd^Db + 



1 



2 cos Tj 



,s s +s s 



Sz)Bi, — sinrjs'^SzFhd 



+ {s-y{d+yEb - smr]SzS-d+Fb + -(sj - Sz)Cb, 



cos(2cj) cos 7] sin(?7 + 2d^'r]) sin(2?7 + Ad'^i]) 



+ [(i sin rj — sui{d^ri) sin(r/ + d^-r])] • [2 cos r] sin(2r7 + 2d^ri) 
4 cos(^ + ^) cos r/ sin(cj + rj + d^rj) sin(— q + 77 + d'^rj) 



+ 



2sin(2^)(d sin 7/ — sin(r7 + d^??) sin(2r/ + d^rj)) 1 
sin(2-e)sin(2+e) 7' 

- — - -^ 2sinrysin(2cj)[d sin r/ — cos?7sin((i^r/)] + sin(2r7)(sin (2(i^r/) — sin q 



2pd(ci) 

Z cos l^^~r^j /j2-2 • -/ \ n ■2/7?\ 

jrp ^ X (a sm r] + sm Cj sin( Cj + 77) — 2 cos r/ sm (a rj) 

sin(^ + 

x(d sin 77 + sin Cj sin(cj — ry) — 2 cos r/ sin {d^rf) 

4sin(2^) cosr7sin(cj — d^rf) sin(— Cj — d^r]){d'^ sin^ r] — sin{d^r]) sin(r/ + d^r])) 

sin(2-e)sin(2+e) 

— - < cos{2ci) cos 7] sinirj — 2d^r]) sin(2r/ — Ad'^r]) 

^Pd{ci) I 2 

+ [d sin ?7 + sin((i^r7) sin(r/ — d^ry)] • [2cosr/sin(2r7 — 2#r/) 

4 cos(^ + ^) cos ry sin(ci + ry — #77) sin(+Ci — 77 + #r/) 



+ 



sin(2+e) 
2 sin(2^) sin(cj — d^rj) sin(ci + 77 — d^rj) sin(ci + ^^77) sin(cj — rj + d^rj) 



sin(2-0sin(2+^) 



sin(cj — ^ — d^ri + ^) 
2sin(2-e)sin(2+e) 



.?? 



,V 



,5r] 



{2d sin 77 — cos rj) sin( ^) — cos(2cj) sin( — h + ^i^ ^7 '^os( 1" '^d^r] — ^) 



'2 '' V w ^2 -' ' ' 2 

sin"^ 77sin(cj — ^ — d^?? + C) sin(ci — 2 ~ d^V + 



sm Ci 



2 cos 77 sin(2 - ^) sin(2 + ^) 
- d^77 + 



2sin(2-^)sin(2+e) 



-{2d sin ?7 — cosr/) sm(— + ^) + cos(2ci) sm(- — ^) — sin?7cos(— — 2(i^77 + ^) 

where i = a,b, and ^ should be changed to ^^ when i = b. As in the spin-1/2 case, one can 
check that the Hamiltonian is hermitician when we choose pure imaginary q. 

3.4 The Bethe ansatz for this model 



To construct the algebraic Bethe ansatz, we define the pseudo- vacuum state |0) as 
st\0) = d+|0) = 0, {i = l, 2, •••, N) 
#|0) = m|0). (43) 

And as before, we write T(u) as 

■ Aiiu) Bi{u) B2{u) 
T{u) = I Ci(n) A2{u) B^{u) I . (44) 

C2{u) C-i{u) A^{u) 

In order to simplify our calculation, we introduce the following transformations 

7 / s - . X sin(2r?) , , , , , 

^^<") = -^''">- sin(2u + 2,) -^''"' C^* 

^3(,) = A^iu-i-'-^^Mu)- . ./""'""'f,"' , M u) (46) 

^ ' ^ ' sm(2u) ^ ' sm(2n + 7?) sm(2u, + 2r?) ^ ' ^ ' 

It is easy to show 

Ci(tx)|0)=0, fi,(n)|0)/0, (i = l, 2, 3), 

Aiium = u;i|0), ^2(n)|0) = w;2|0), ^3(^)|0) = w^, 

with 

u^i = sin^ [u + rf) s\v? {u + 2?]) 

X IT IT sin(n + rcj + m?]) sin(ti + rcj + r/ + m?]), 

i=a,6r=itl 

sin(2^)sin(^-K-|^) . ^^ ■ 2N r ^ ^ 
sm(2n + z?7j sm(^ + u + ^ j 
X IT IT [sin(ii + rcj) sin(n + rci + rj) + d sin rj — 2 sin (mr/) cos ry, 

i=a,b r=±l 

sin(-^ + u+ 2)sin(-^ + u+ |^)sin(7?-2u) . ^j^ . ^j^ 

W3 = — r-rz n , ■ / — 7 nt ■ / ^T^ sm u sm (u - rj) 

sm(^ + u+ 2)sm(-^-ti+ 2)sm(r/ + 2ti) ^ '^ 

X IT IT sin (n + rCj — m?]) sin (u + rcj + r/ — mry). 
i=a,6r=±l 

From the reflection equation (5), we can find the following commutation relation 

, / N^ / M V sin(u — f — 2r7) sinfii + w) ^ , . . , ^, , 

Aiiu)Biiv)0) = — ^ ^i . ) {Biiv)Aiiu)0) 

^ ^ -^^ ^1 / sm(n + t; + 2ry) sm(u - w) ^ ^ ^^ ^i / 

sin(2?;) sin(2r7) ^ ,,.,-,, , 
+ — 7 \ / -Bi(u)Ai(v)0) 

''''^^''^ Bi{u)A2{vm (47) 



sin(u + V + 2ri) 
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A2iu)Biiv)\0) 

sm{u — V — 2rj) sin(u — v + rf) sin(ti + v) sin(ti + v + 3r]) 
sin(n — V — rf) sin(u + v + rf) sin(u — v) sin(u + v + 2rf) 



Bi{v)A2{u)\Q) 



+ —( \ - V9 ^9 M n)Ai{v)\^)- . , ^'' M u)A2{vm 

sm(ti — V — r}) sin(2f + 277) sva.[u + f + r/j 

I sin(27?)[sin(2T?) sin(if - v + r]) - sin(n + i; + t?) sm(2u + 2t?)] , . t . .,„, , .„. 

+ ; ^ -, ^ -, ^ Bi(u)A2\v)\yj) (48) 

sin(27/) sin(2t;)[sin(2n + 3r/) sm(ti — f — 2r/) — smrysin(n + w + 2r/)] 

+ ^-7 T— 7 — \ , o A ■ ^o — I n \ ■ /o — 7^^\ Oi[u)Ai{v)\\)), 

sin(n — V — rj) sin[u + v + 2r]) sin[2u + 2ri) sin(2t; + 27]) 



A3iu)Bi{v)\0) 

sm{u — V + r]) sin(n + v + 3r/) 
sin(n — V — rf) sin(u + w + 77) 
s\n.{2rf) sin(2f ) sin(2u + 27^) 



Bi{v)Az{u 



sin(2ii) sin(ti — v — rf) 



B3{u)A2{v)\Q) 



sin(2r?)sin(2.)sin(2n + 2r,) ^ (^)^ (^)|o^ 
sin(2n) sin(2v + 2r/) sm(ii + t; + r/) ^ ^ ^v ;i / 

sin^(277) sin(2n + 2r?) ^ . x 7 , x, > 

sm(2n) sin(2ti + rj) sin(n — v — rj) 

sin2(2,) sin(2.) sin(2n + 2,) ^,(,)^,(,)|o), (49) 



sin(2n) sin(2ti + r/) sin(2t; + 2r]) sm{u + v + 2rf) 
From the reflection equation (5), we can constructed a two-particle excited state as follows, 

\V^,V2) = {b^{vi)Bi{v2) + -. ''''^^'^^ , B2{V1)A2{V2) 

I sm vi + sm V2 + ri 

- ■ /^^(^^)^^"(;^ + ^^-^) n,^,,)A{v2)]\^) (50) 

sin(T;i -V2-rj) sin(z;i +V2 + rj) J 

which is symmetric in ^1,^2 up to a whole factor. Applying the transfer matrix on the two- 
particle excited state, we have found a lot of "unwanted terms" . They must cancel each other 
to ensure the above state to be eigenstate. However, we can't check it directly because the 
calculation is much more complicated than we expected. Here we assume that they vanish. 
To get the Bethe ansatz equations, we have to use the functional Bethe ansatz method which 
was first proposed for the Ising model|2^. Tarasov et al argued that this method can be 
generalized for n-particle excited state[29-31]. Then from equation (8), the eigenvalue of the 
transfer matrix for n-particle excited states is as follows 

i{u)\vi,V2,--- ,Vn) = {kfAl{u) +k^A2{u) + k^A3{u))\vi,V2,---,Vn) 
= {w^Al{u) +w:^A2iu) +W^A3{u))\vi,V2,---,Vn) 

, -J^ sm{u - Vi - 2ri) sm{u + Vi) 

= Wlwi[[^- ■ , ^ X ■ . TVi,V2,---,Vn) 

^_^ sin(u + Vi + 2rj) sin[u — Vi) 
, -A- sin(ii — Vi — 2r/) sin('u — Vi + r]) sin(u -|- Vi) sin(u + Vi + 3r]) 

+ WJW2[[^^ :——, ■ ■ :——, ^— ^ ■ ——Vl,V2,---,Vn) 

V: sm(n — Vi — r]) sin[u + Vi + r]) sin[u — Vi) sin(n + Vi + 2rj) 

, -^ sm.{u - Vi + rf) sm{u + Vi + ?,r]) 

+ wJw^U^— — — — ■ ■ — Avi,V2,--- ,Vn) (51) 

^J- sm(u — Vi — r]) sin[u + Vi + r]) 
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where 



Wi [u 



Wn [U] 



w^{u) 



sin(2n + 3r/) sm{u — ^"' 



sin(2ii + r/) sm(u — ^+ + -^) 

sin(^+ +u + ^) sm{u - ^+ - 2) sin(2u + 2r/) 
sin(-C+ + u + ^) sin(-n + e+ - i) sm(2w) ^ 
sin(e+ + u+^) sm{-u - g+ - ^) 

sin(-^+ + n + |i) sin(-n + C+ - ^^ ' 



and the free parameters fj,i=l,2,- • -,71, obey the Bethe ansatz equation 

sin{^+ - Vj + '^) sin{^ + Vj - ^) sin'^'^ {vj + ri) -^ -p-j- sm{v j + rck + mr]) 

sin(^+ + Vj + ^) sin(^ - Vj - D sin2^(i;j - r/) ^^^^^,^±1 ^^"^(^J + '^'^'^ ~ "^'^) 
= TTsinfa-.. + r?)sin(.,+.. + r?) , = i, 2, • • • , n. (52) 



«7^i 



sin(7;j — Vi — 7]) sm{vj + Vi — r])'' 



This Bethe ansatz equation can also be derived by the means of the fusion method, we 
have checked that they agree with each other 1 32, R^. With this method of Bethe ansatz, we 



cannot get complete eigenstates. However, it is a powerful tool to obtain the eigenvalues and 
Bethe ansatz equation for models which cannot be obtained with fusion method [pS|]. 

4. Discussion 

In this paper, we studied the integrability of the spin-1/2 and spin-1 XXZ open Heisen- 
berg chains with boundary impurities. These models are relevant to the Kondo problem in 
a Luttinger liquid. By using the algebraic Bethe ansatz and its extension, we have obtained 
the eigenvalues of the Hamiltonians and the Bethe ansatz equations. When we let d = 
in this paper, one can easily check that the Li{u) {i = a, b) in formulas (10) and(40) will 
be identity so the present Hamiltonians and Bethe ansatz equations can be reduced to the 
usual ones respectively. This procedure can be generalized to the general Heisenberg chain. 
It is worthy to point out that the Bethe ansatz equations and eigenvalues of transfer matri- 
ces for the spin-1/2 and spin-1 XXX chains coupled with arbitrary spin impurities can be 
obtained by rescaling the spectral parameters Vj by Vj x r] and taking the limit r/ — > in 
Bethe ansatz equations and the eigenvalues of the transfer matrices. With similar methods 



of Ref . [p4|, |35| , |36| , the results of the present paper can also be used to calculate the boundary 
susceptibility, the contribution of the impurities to the specific heat and Kondo temperature, 
which can describe the effect of impurities to the system. We will study them in another 
paper! 
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